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Abstract 

We study the properties of the G + in nuclear matter and + hypernuclei within the quark 
mean-field (QMF) model, which has been successfully used for the description of ordinary nuclei 
and A hypernuclei. With the assumption that the non-strange mesons couple only to the u and 
d quarks inside baryons, a sizable attractive potential of the + in nuclear matter is achieved as 
a consequence of the cancellation between the attractive scalar potential and the repulsive vector 
potential. We investigate the + single-particle energies in light, medium, and heavy nuclei. More 
bound states are obtained in G + hypernuclei in comparison with those in A hypernuclei. 
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I. INTRODUCTION 



The signature of the pentaquark was found in the experiment with GeV photons on 12 C at 
LEPS of SPring8 Q, ^. The jn — > K + K~n reaction was analyzed and a narrow resonance 

state was identified with the mass of 1.54 GeV and the width smaller than 25 MeV. Such a 

n 

narrow state was predicted by Diakonov et al. [3( as a pentaquark state at nearly the same 
energy. There were many experimental works published after the original study by using 
various reactions |4j. Recently, many further experiments, particularly at very high energy, 
report negative results except for a few positive results |4j. At this moment, the spin and 
parity of the state are not known. Even the width is not precisely identified experimentally. 

There are many theoretical works [4]. It seems those models, which emphasize the im- 
portance of the chiral symmetry, suggest the positive parity (J p = l/2 + ) as the original 
predication of Diakonov et al. [3|. On the other hand, the naive constituent quark models 
tend to provide the negative parity (J p = l/2~). The QCD originated works as the lattice 
QCD and the QCD sum rule approaches seem to provide the negative parity, although these 
methods face the difficulty of handling the nearby K + n threshold and the small pion mass. 
Hence, we are at this moment not sure about the presence of the pentaquark state, the 
spin-parity and even the width of this state 0]. 

We should make further efforts to get more information for the pentaquark state. There 
are several interesting theoretical works to study the pentaquark in nuclear medium 0, 0, Q, 
0, |^ . In Ref . j3| , the authors study the interaction of the pentaquark with nucleons using the 
meson exchange model. They have obtained a large binding potential with the width to be 
small to take the pentaquark nucleus seriously. This work then motivated a theoretical work 
to calculate the (K + ,tt + ) reaction for the formation of the pentaquark nuclei Q|. There is 
a good possibility to see + hypernucleus in such a reaction. It is then very interesting to 
study the + hypernucleus from different view point. 
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We studied nuclei and hypernuclei in terms of the quark mean-field (QMF) model 
We assume in the QMF model that up and down quarks in the baryons interact with a, 
u, and p mesons (mesons made of up and down quarks) in nuclear matter and in nuclei. 
Hence, the interaction of the + with other nucleons is obtained without introducing new 
parameters in the QMF model. The only assumption needed here is the isospin of the G + 
to be / = 0. The interactions of quarks with the meson fields change the properties of the 
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baryons in hadronic matter and hence the hadronic matter property is obtained from these 
interactions self-consistently. 

In the next section, we would like to present the formulation of the G + hypernuclei in 
the relativistic quark mean-field model. In section HTT1 we provide the numerical results for 
various + hypernuclei. Section |TV| will be devoted to the conclusion and discussions. 



II. QUARK MEAN-FIELD MODEL FOR THE 0+ IN MEDIUM 

We start with the description of the G + in nuclear medium, which is comparable to the 
treatment of A hypernuclei in our previous work Q]. The G + is known to be an exotic 
baryon containing five quarks (uudds) with strangeness S — +1, isospin 1 = 0, and charge 
Q = +1. However, its parity is still ambiguous experimentally. In this work, we use the 
constituent quark model to describe the B + and nucleons in medium, while the exchanged 
mesons couple directly to the quarks inside baryons. The constituent quarks satisfy the 
Dirac equation: 

-m q - X c- 9l° - fl£w7° - <?>r 3 7°] 9 (r) = 0, (1) 

where \c is the confining potential taken as the quadratic form, \c = \kr 2 . Because the 
Lorentz structure of the confining interaction is not well known yet, we here use a scalar 
confinement. In our previous studies for finite nuclei and A hypernuclei [3, 12], both 



scalar and scalar- vector confinements have been taken into account, but it is known that 
antiquarks cannot be confined by a scalar-vector confinement when the vector part is equal 
or larger than the scalar part. Therefore we just adopt the scalar confinement in the present 
calculation of G + . We assume that the non-strange mesons, a, u, and p, couple exclusively 
to the up and down quarks and not to the strange quark (or antiquark) according to the 
OZI rule, hence the strange antiquark inside B + satisfies the Dirac equation in which the 
terms coupled to the non-strange mesons vanish (g^ s = g^j s = g s p ,s = 0). We note that the 
Dirac equation for the strange antiquark s in + is identical to the equation for the strange 
quark s in A when a scalar confining potential is adopted in Eq. We follow Ref. ll( to 
take into account the spin correlations and remove the spurious center of mass motion, and 
then obtain the effective mass of the + in medium as 



Me = J(*e q + e- s + £f pin ) 2 - (4(rf> + (pg», (2) 



while the effective mass of the nucleoli is expressed as 



M* N = ^(3e q + E» in y-3(pt). (3) 

Here, the subscript q denotes the u or d quark. The energies (e q and e s ) and momenta {(p 2 ) 
and (pf)) can be obtained by solving the corresponding equations. The spin correlations 
(i?® in and E^ in ), which are taken as parameters in the QMF model, are determined by 
fitting the baryon masses in free space (Me = 1540 MeV and Mn = 939 MeV). The 
effective masses of baryons are influenced by the a mean-field in nuclear medium, which 
provides a scalar potential to the u and d quarks in baryons and as a consequence reduces 
the constituent quark mass to m* = m q + g q a o (q = u,d). However, the u and p mean fields 
could not cause any change in the baryon masses, and they appear merely as the energy 
shift. Therefore, we obtain the effective masses of baryons as functions of the a mean-field, 
M e (<r) and M* N (a), expressed by Eqs. and ©. 

We now derive a self-consistent treatment for a G + hypernucleus, which contains a pen- 
taquark + and many nucleons. The baryons in this system interact through the exchange of 
a, u, and p mesons, just like in the case of a single A hypernucleus. The effective Lagrangian 
at the hadron level within the mean-field approximation can be written as 



£ = ip 



iry^ - M* N (a) - - g P pr 3 l° ~ efi^pUy 



i> (4) 



- Ml (a) - - eA^] 



-\(V°) 2 ~ \rn 2 y - ^g 3 a A + ^(v^) 2 + \™l^ 2 + \c 3 u A 

where ip and ip® are the Dirac spinors for the nucleon and + . The mean-field values are 
denoted by a, u, and p, respectively, while m a , m w , and m p are the masses of these mesons. 
The electromagnetic field, denoted by A, couples both to the proton and B + , which carry a 
positive charge. The influences of the a mean-field on baryons are contained in the effective 
masses M^(cx) and M (cr). The uo meson, which couples directly to the u and d quarks inside 
baryons, provides an interaction at the hadron level with the coupling constant = 
where Tij is the number of the u and d quarks in the baryon i. When % = + (z = N), we 
get g® = ig^j (g^ — — 3g%). The p meson doesn't couple to the + which is an isoscalar 
particle, but it couples to the nucleon with the coupling constant gf? — g p — g% as given 



in Ref. [11 1. In the QMF model, the basic parameters are the quark-meson couplings (g%, 
g%, and the nonlinear self-coupling constants (#3 and C3), and the mass of the a meson 
{m a ), which have been determined by fitting the properties of nuclear matter and finite 
nuclei in Ref. Therefore, no more adjustable parameters exist when it is extended to 
the calculations for B + or A hypernuclei. From the Lagrangian given in Eq. (@J), we obtain 
the following Euler-Lagrange equations: 

1 + 7*) 



i TM 0" - M* N (a) - - g p pr^ - e- 



-A7 



^ = 0, 



i lv & - M* (a) - g®^ - eAj } ^ e 

P. 



0. 



-A + vnf a ) a 

,2' 



dM* N (a) 



da 



9M* e (a) e 

rs 



A + mf) p = g p p 3 , 



da 
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-AA 



(pp + Pv 



(5) 

(6) 

(7) 
(8) 
(9) 
(10) 



where p s (pf), p v (pf), P3, and p p are the scalar, vector, third component of isovector, and 
proton densities, respectively. The above coupled equations are solved self-consistently with 
the effective masses M^(a) and M@(cr) obtained at the quark level. 

The application of the QMF model to the description of A hypernuclei has been reported 



in our previous work [12j. Without adjusting any parameters, the properties of A hyper- 
nuclei could be described reasonably well. The A single-particle energies obtained in the 
QMF model are slightly underestimated in comparison with the experimental values. In the 
present work, we investigate the + hypernuclei within the QMF model, and then provide 
the predictions for + single-particle energies in light, medium, and heavy nuclei. 



III. NUMERICAL RESULTS 



In this section, we will present the results calculated in the QMF model for the + in 
nuclear medium. We need first to specify the parameters used in the present calculation. For 
the parameters at the quark level, we take the constituent quark masses m u = rrid = 313 MeV 
and m s = m s = 490 MeV, and the strength of the confining potential k = 700 MeV/fm 2 , 
as those used in the case of A hypernuclei The meson masses are taken as m a = 

470 MeV, = 783 MeV, and m p = 770 MeV. The quark-meson couplings, g q a = 3.14, 
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g% = 4.20 and g q p = 4.3, and the nonlinear self-coupling constants, = 50.7 and C3 = 53.6, 
have been determined by fitting the following equilibrium properties of nuclear matter: 
equilibrium density po — 0.145 fm -3 ; binding energy E/A = —16.3 MeV; incompressibility 
k = 280 MeV; effective mass = 0.63 M^; symmetry energy a sym = 35 MeV. With the 
above parameters, the QMF model has provided quite satisfactory results for finite nuclei 
and A hypernuclei Now, we take the same parameters as above to perform the 

calculations for the + in nuclear medium. 

It is very interesting to see how the G + properties change in nuclear matter. In Fig. 1, 
we present the variation of the effective mass of the B + , in comparison with those for the 
nucleon and A, as a function of the quark mass correction due to the presence of the a 
mean-field, 5m q = m q — m* = —g q a o (q = u,d). It is obvious that the reduction of Mq is 
larger than those of Mjy and M£, since there are four u and d quarks in the + which are 
influenced by the a mean-field, but only three or two quarks influenced in the nucleon or A. 
We note that the dependence of the effective masses on the a mean-field must be calculated 
self-consistently within the quark model, therefore the effective masses of baryons obtained 
in the QMF model are not simply linear functions of a, as given in the relativistic mean- 



field (RMF) models 13J, [14J. We show in Fig. 2 the ratios of the baryon radii in medium 



to those in free space as functions of the nuclear matter density. It is found that all of the 
baryons, + , A, and N (nucleon), show significant increase of the radius in nuclear medium, 
which is consistent with the famous EMC effect. The nucleon gains larger increase of the 
radius than the G + and A, because all three quarks in N are involved in the interactions 
with mesons while only four-fifths in + and two-thirds in A involved. The baryon radii 
increase by about 4% ~ 6% at normal matter density in the present model. In Fig. 3, 
we plot the scalar and vector potentials of the B + in medium as functions of the nuclear 
matter density, and the results of the nucleon (N) and A are also shown for comparison. At 
p = p = 0.145 fur 3 , we get Uf = -420 MeV and U§ = 370 MeV, therefore an residual 
attractive 9 + potential, U e = Uf + Uy = —50 MeV, is predicted in the present model. 
For A at p = po, we have got the potential U A = Ug + Uy = —24 MeV, which is about 



half of the (~) + potential. In Ref. Q], a large attractive B potential ranging from —GO to 
— 120 MeV at normal matter density was predicted by calculating the G + selfenergy tied 
to the KN decay of the + and the two meson baryon decay channels of the + partners 
in an antidecuplet of baryons. A QCD sum rule calculation ,8| predicted an attractive + 
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potential of about —40 ~ —90 MeV. Therefore, the prediction of the + potential in the 
QMF model is compatible with the results obtained by other groups [?], [| . 

We now present the results of self-consistent calculations for the G + hypernuclei con- 
taining a + bound in nuclei. In Fig. 4, we show the predicted + single-particle energies 
in q 7 0, gCa, and Q^Pb. It is found that there are many bound states of the + in light, 
medium, and heavy nuclei. We can see that the separation between the deep G + energy 
levels is comparable with the + width. The theoretical prediction of the free B + width 
in the chiral soliton model is less than 15 MeV Q], and the + width in medium might be 
reduced to about one third or less of the free width due to Pauli blocking and binding, as 
pointed out in Ref. It is clear that large separation of the deep + levels is helpful to 
observe distinct peaks in the experiment to produce G + bound states ^(|. It is very inter- 
esting to compare the results of G + hypernuclei to those of A hypernuclei. We can see that 
the single-particle energy of lsi/2 @ + in Q 09 Pb is smaller than the one in @Ca, which is con- 
trary to the results of A hypernuclei. This is mainly because @ + carrying a positive charge 
behaves like a proton in nuclei, while A is a charge neutral particle more like a neutron. 
The Coulomb energies of G + in heavy nuclei are much larger than those in light nuclei, that 
+ single-particle energies in heavy nuclei can be reduced by the large positive Coulomb 
energy. We plot in Fig. 5 the scalar, vector, and Coulomb potentials in @Ca with a @ + 
bound at lsi/2 state. At the center of the G + hypernucleus, the attractive scalar potential 
(Ug ~ —559 MeV) is partly cancelled by the repulsive vector potential (Uy ~ 515 MeV) 
and Coulomb potential (Uc ~ 11 MeV), and as a consequence, an attractive 6 + potential 
(U ~ —33 MeV) remains in this case. 

It is interesting and important to compare the results of the QMF model with those ob- 
tained by other groups and discuss the origin of the differences in the B + single-particle 
energies. More bound states of the G + in nuclei were predicted in some recent works j?l E| ■ 
In Ref. j?| , the authors obtained more deeply bound states by solving the Schrodinger equa- 
tion with two potentials: V(r) = -60p(r)/p (MeV) and V{r) = -120p(r)/p (MeV). It 
is obvious that larger 6 + single-particle energies should be achieved due to the deeper + 
potentials used in their calculations. In Ref. 9j, the + hypernuclei were studied by using 
the RMF model with the 9 + couplings, gf = ~g a and g® = ^g^, which are the quark model 
predictions close to those given in the quark-meson coupling (QMC) model. It may be 
helpful to look at the case of A hypernuclei, because of the similarity of the B + hypernuclei 
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and A hypernuclei. It is well known that the properties of A hypernuclei are very sensitive 
to the effective coupling constants at the hadron level, especially the two relative couplings 
R% = g^/g a and = g^/g^ [15]. The quark model prediction, = R^ = 2/3, usually 
gives large overbinding of A single-particle energies in comparison with the experimental val- 
ues. For example, by taking TM1 parameter set in the RMF model with R% = R^ = 2/3, 
the lsi/2 A single-particle energy in ^Ca is calculated to be —35.2 MeV, while in the case 
of the NL-SH parameter set used, it is —36.8 MeV. We note that the experimental value 
of the lsi/2 A binding energy in ^°Ca is about —18.7 MeV. Usually some small deviations 
from R% = R^ = 2/3 should be taken into account in the studies of A hypernuclei in order 
to reproduce the experimental values [3]. It has been pointed out that about 10% increase 
from R~ = 2/3 seemed to be needed to improve the overestimated A single-particle energies 
in the QMC model, since the R% calculated at the quark level in the QMC model was very 



close to 2/3 16, 17]. On the other hand, the QMF model yielded slightly underestimated A 
single-particle energies due to smaller value than 2/3 of the R% obtained in the QMF model, 
which required about 3% reduction from R^ = 2/3 in order to reproduce the experimental 
values. It should be due to the same reason that the QMF model predicts relatively small @ + 
single-particle energies in comparison with those obtained in the RMF model with the quark 
model value of the 9 + coupling, gf / g a = g® /g^ = 4/3. In the QMF model, g® = 4g^ = -g^ 
is adopted according to the assumption that non-strange mesons couple only to the u and 
d quarks inside baryons, but gf = has to be worked out self-consistently at the quark 
level. We notice that the effective a— baryon couplings in the QMF model are not simply 
constants as in the RMF model. Therefore, the relative coupling obtained in the QMF 
model, R® = —j^- / > depends on the a mean-field in medium, and we get smaller 

values than 4/3 for the R® in the QMF model. This is considered as the dominant origin of 
the differences between the predictions given by the QMF model and those obtained in the 
RMF model Q. 

In order to examine the sensitivity of + single-particle energies to the + couplings, 
we perform the same calculation with a reduced coupling 0.97 x g®. In this case, the lsi/2 
@ + single-particle energy in ^Ca changes to be —38.7 MeV from the value —28.1 MeV 
as shown in Fig. 4. Hence, 3% reduction of g® can lead to ~ 10.6 MeV decrease of the 
lsj/2 + single-particle energy in ^Ca, and it provides ~ 10 MeV more attraction for the 
@ + potential in nuclear matter at normal matter density. In the present work, we have 
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performed self-consistent calculations for both 6 + and nucleons in + hypernuclei. By 
comparing the results of a normal nucleus containing A nucleons with those of the + 
hypernucleus containing A nucleons and one 6 + , the effects of the 6 + on nucleons can be 
found. The existence of the + does enlarge the scalar and vector potentials of nucleons, 
but the enlargements are mostly cancelled with each other. 

IV. CONCLUSIONS 

We have studied the properties of the 6 + in nuclear matter and performed self-consistent 
calculations for 6 + hypernuclei in the framework of the QMF model, which has been suc- 
cessfully applied to the descriptions of stable nuclei and A hypernuclei. In the present work, 
we have used the constituent quark model to describe the + and nucleons, which natu- 
rally allows the direct coupling of exchanged mesons with quarks inside baryons. With the 
parameters determined by the equilibrium properties of nuclear matter and the assumption 
that non-strange mesons couple only to the u and d quarks, an attractive + potential has 
been achieved as a consequence of the cancellation between the attractive scalar potential 
and the repulsive vector potential in nuclear matter. The QMF model provides an attractive 
@ + potential of ~ —50 MeV, and yields a reduction of effective + mass of ~ 420 MeV as 
well as ~ 5% increase of the G + radius at normal matter density. 

We have investigated the 9 + single-particle energies in light, medium, and heavy nuclei. 
It is found that there are many bound states in + hypernuclei, and the 6 + seems to be 
more deeply bound than the A in hypernuclei. We have discussed the sensitivity of + 
single-particle energies to the + couplings. The comparison between the results calculated 
in the QMF model and those obtained by other groups has been done, and the origin of the 
differences has also been discussed by referring to the case of A hypernuclei. The results 
obtained in the QMF model provide encouragement for the experimental work to produce 
6 + hypernuclei. 

We comment here a possible method to produce 6 + hypernuclei. The gamma induced 
reaction on nucleus for the formation of the 6 + leads to the 6 + in a scattering state due to 
the large momentum transfer. On the other hand, the K + induced reaction with a nucleus 
could produce the + with a reasonable momentum of order 500 MeV, when an outgoing 
pion is observed at the forward angle. This momentum transfer would lead to high angular 
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momentum transfer of order L ~ 10, when the + is produced at the nuclear surface with a 
medium heavy nuclei. Hence, it is highly plausible to observe a sharp peak with high angular 
momentum; high spin B + - high spin nucleon hole state, popping up in the beginning of the 
continuum region as the case of the (p,7T~) reaction [l^ . 

Acknowledgments 

This work was supported in part by the National Natural Science Foundation of China 
(No. 10135030) and the Specialized Research Fund for the Doctoral Program of Higher 
Education (No. 20040055010). 



10 



[1] T. Nakano et al. (LEPS collaboration), Nucl. Phys. A721, 112c (2003). 

[2] T. Nakano et al. (LEPS collaboration), Phys. Rev. Lett. 91, 012002 (2003). 

[3] D. Diakonov, V. Petrov, and M. Polyakov, Z. Phys. A 359, 305 (1997). 

[4] http://www.rcnp.osaka-u.ac.jp/~penta04/ 

[5] G.A. Miller, Phys. Rev. C 70, 022202 (2004). 

[6] H.C. Kim, C.H. Lee, and H.J. Lee, J. Kor. Phys. Soc. 46, 393 (2005). 

[7] D. Cabrera, Q.B. Li, W. Magas, E. Oset, and M.J. Vicente Vacas, Phys. Lett. B 608, 231 
(2005). 

[8] F.S. Navarra, M. Nielsen, and K. Tsushima, Phys. Lett. B 606, 335 (2005). 

[9] X.H. Zhong, Y.H. Tan, L. Li, and P.Z. Ning, Phys. Rev. C 71, 015206 (2005). 

[10] H. Nagahiro, S. Hirenzaki, E. Oset, and M.J. Vicente Vacas, | nucl-th/ 0408002 

[11] H. Shen and H. Toki, Phys. Rev. C 61, 045205 (2000). 

[12] H. Shen and H. Toki, Nucl. Phys. A707, 469 (2002). 

[13] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

[14] Y. Sugahara and H. Toki, Prog. Theor. Phys. 92, 803 (1994). 

[15] CM. Keil, F. Hofmann, and H. Lenske, Phys. Rev. C 61, 064309 (2000). 

[16] K. Tsushima, K. Saito, and A.W. Thomas, Phys. Lett. B 411, 9 (1997); (E) Phys. Lett. B 
421, 413 (1998). 

[17] K. Tsushima, K. Saito, J. Haidenbauer, and A.W. Thomas, Nucl. Phys. A630, 691 (1998). 

[18] A.B. Brown, O. Scholten, and H. Toki, Phys. Rev. Lett. 51, 1952 (1983). 



11 




FIG. 1: The variations of the baryon effective masses, SM* = M* — Mi (i = N, A, + ), as functions 
of the quark mass correction, Sm q = m q — m* = —g%a. 
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FIG. 2: The ratios of the baryon radii in nuclear matter Ri to that in free space as functions 
of the nuclear matter density p. 
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FIG. 3: The scalar and vector potentials, Us and Uy, as functions of the nuclear matter density p. 
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FIG. 4: The + single-particle energies in q O, @Ca, and @ 09 Pb. 
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